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Abstract 

We describe an imbalanced superfluid Fermi gas in three dimensions within the path-integral 
framework. To allow for the formation of the Fulde-Ferell-Larkin-Ovchinnikov-state (FFLO-state), 
a suitable form of the saddle-point is chosen, in which the pairs have a finite centre-of-mass mo- 
mentum. To test the correctness of this path-integral description, the zero-temperature phase 
diagram for an imbalanced Fermi gas in three dimensions is calculated, and compared to recent 
theoretical results. Subsequently, we investigate two models that describe the effect of imposing a 
one-dimensional optical potential on the 3D imbalanced Fermi gas. We show that this ID optical 
potential can greatly enlarge the stability region of the FFLO-state, relative to the case of the 
3D Fermi gas without ID periodic modulation. Furthermore it is show that there exists a direct 
connection between the centre-of-mass momentum of the FFLO-pairs and the wavevector of the 
optical potential. We propose that this concept can be used experimentally to resonantly enhance 
the stability region of the FFLO-state. 
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I. INTRODUCTION 



The study of superconducting and superfluid systems has recently attracted wide at- 
tention, among other things because of the realization of ultracold Fermi gasses in optical 
lattices jl-4]. These systems can be considered as quantum simulators that can be used 
for probing fundamental problems in condensed- matter physics [5(, for instance the search 
for exotic new phases in strongly magnetized superconductors. Ultracold Fermi gasses offer 
important advantages over conventional superconductors, mainly because of their extensive 
tunability. In a superconductor, the number of spin-up and spin-down electrons is equal and 
the interaction strength is fixed. In ultracold Fermi gasses, one can not only tune the inter- 
action strength by use of Feshbach resonances 6-8j], but also the population imbalance can 
be freely adapted. This experimental freedom has led to the study of a variety of new phe- 
nomena in imbalanced ultracold Fermi gasses (9 13 ] . One fundamental question, that is still 
not settled, concerns the nature of the ground state of an imbalanced Fermi gas. When pop- 
ulation imbalance between the spin-up and spin-down components is introduced into these 
systems, complete pairing is no longer possible. Clogston and Chandrasekhar suggested that 
above a critical imbalance, the superfluid system would undergo a transition into the norma 
state [14], [l5|. This effect has been observed experimentally by the MIT 9] and Rice [l0| 
groups. However, their observations were not in exact agreement, and there still exists some 
controversy [16J about the exact nature of the phases of the superfluid system at high levels 
of imbalance. In 1964 Fulde and Ferell and independently Larkin and Ovchinnikov [lsl ] 
proposed that a superfluid system can accommodate population imbalance, by making a 
transition into a state with a finite center-of-mass momentum (and thus a spatially mod- 
ulated order parameter). This state is the so-called Fulde- Ferell-Larkin-Ovchinnikov-state 
FFLO-state). Recently, there has been an ongoing theoretical search for this exotic state 
JiJ. In ID, the FFLO-state was predicted to exist in superconducting systems 20]. This 



work was confirmed numerically 21[ and elaborated further 



irough theoretical studies of 



the ground state and the phase diagram of a ID Fermi gas |22M24I] . Furthermore it has been 
shown that the ID analogue of the FFLO-state is stable in a large section of the BCS-BEC 
crossover phase diagram 25), compared to the case of a 3D Fermi gas. Although this ID 
FFLO-state has not yet been observed directly, a recent paper reports the experimental 
observation of density profiles that agree quantitatively with theoretical predictions at low 
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temperature 



261 ] . In three dimensions however, the experimental observation of the FFLO- 



state has so far remained elusive. One of the main reasons for this is that the FFLO-state 
in three dimensions only occurs in a tiny section of the BCS-BEC crossover phase diagram 



n 



This then begs the question, is there a way to stabilize the FFLO-state in a 3D 
Fermi gas? The purpose of this paper is twofold: first we develop a path-integral description 
for a superfluid Fermi gas which can accommodate the FFLO-state, and second we propose 
a method to stabilize the FFLO-state through an optical potential. In two recent papers it 
was suggested to stabilize the FFLO-state by the use of a 3D optical lattice 29|, |39|]. In this 
paper we investigate the stabilizing effect of a ID optical potential in order to investigate 
the interplay between the wavevector of the FFLO-state and the wavevector of the laser 
which creates the optical potential. In the present work, the ID optical potential provides a 
periodic modulation in one direction. We emphasize that we do not look at the FFLO-state 
in a one-dimensional gas [30] , but in a 3D gas with a superimposed one-dimensional periodic 
potential. In section [TT] we describe the FFLO state in an imbalanced Fermi gas in 3D within 
the path-integral framework. As a test for the correctness of this description, we calculate 
the zero temperature pha se diagram for this system and compare our findings with recent 



theoretical results 



27 



28|. In section HT] we investigate two models to account for the effect 



of a one-dimensional optical potential. We show that the presence of such a potential leads 
to a substantial increase of the stability region of the FFLO-state. Finally in section HVl we 
draw conclusions. 



II. PATH INTEGRAL DESCRIPTION 

The partition sum of an imbalanced Fermi gas in 3D can be written as a path integral 
over the fermionic fields ^k,w n ,o- and i/)k tU „,o-' 



Z = J VipVip exp I - y^y^ k,u) n ,<7 (-iu n + k 2 - ipk^ ni<7 

V k,n a 
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/3L 3 

k,n k',ri' q,m 

(1) 

Here k is the wavevector, u n are the fermionic Matsubara frequencies and a =f, -I denote the 
two different hyperfine states. Furthermore, /3 is the inverse temperature given by l/Zc^T, L 
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is the lateral size of the system, the chemical potential of a particle with spin a is denoted by 
ficr and g is the renormalized interaction strength. We use units such that h = 2m = Ep = 1. 
This partition sum (JTJ) can be made more tractable by introducing the Hubbard- Stratonovic 
transformation, which decouples the fourth-order interaction term into second order terms 
by introducing two auxiliary bosonic fields q ,o m and q ,Q m , interpreted as the pair fields. 
As a first approximation, only the saddle-point is taken into account in the path integral 
over the bosonic fields. To describe the FFLO-state, we propose to use a saddle point at 
which the atomic pairs have a finite wavevector Q: 



- , q,n r , 



(S qiQ ^,oV^A . (2) 



By using this particular form of the saddle-point, we choose to describe the FF-state, which 
has an order parameter given by a plane wave ~ e tClr . It is also possible to describe the 
LO-state, which is a superposition of two plane waves with wavevector Q and — Q. In this 
paper, we will not consider the LO-state. For the remainder of the article, the FF-state is 
referred to as the FFLO-state. When Q is set equal to zero in (J2j), the description of the 



normal superfluid is recovered |3jJ. In expression ([2]), the prefactor a/ (3L 3 ensures that A 
has units of energy. Using (T5]) the fermionic fields can be integrated out in expression (pQ) 
for the partition function, leading to an effective action S sp , through 

Z sp = exp ( ln [- det (-Gj^)] + ^ | A| 2 ) = exp (S sp ) , (3) 

V k,n 9 J 

with Gj^ the inverse Nambu propagator which is given by 

-iu n + (Q/2 + k) 2 - m -A 



-A* -iu n - (Q/2 - k) 2 + m 



(4) 



Expression (J3j) can be simplified further by performing the sum over the Matsubara frequen- 
cies. Also, it is useful to express the results as a function of the total chemical potential 
\i = (/Xf + jt/j.) /2 and the imbalance chemical potential ( = (jj^ — /ij,) /2. Furthermore, the 
interaction between particles is modeled with a two-body contact potential V (r) = gS (r). 



The renormalized interaction strength g can then be written as follows 

1 1 ^1 



32|: 



g 8n(k F a s ) ^ 2k 2 ' 



(5) 
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where a s is the 3D s-wave scattering length. As a final step the continuum limit is taken, 
and a thermodynamic potential is associated with the effective action S sp = f3Q sp . This then 
results in 

|2 ' 



n 



up 

U 



^ f i In [2 cosh (/3C Qlk ) + 2 cosh (/3£ k )] - £ Q ,k - ^ 



A 



8ir (kpa s 



(6) 



where the following notations were introduced 

k 2 - 



^Q,k 

/; k <^Q,k 
C Q ,k = C + Q ■ k 
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Cc 2 



A| 



(7) 



The resulting form of the thermodynamic potential has a similar form as the original 
result for the homogeneous 3D Fermi gas derived by Iskin and Sa de Melo 33J| and coincides 
with it for Q — » 0. In order to test the correctness of ([6]), the thermodynamic potential 
is used to calculate the zero temperature phase diagram of an imbalanced Fermi gas in 
3D. This can be done for a fixed number of particles or for fixed chemical potentials. To 
transform between these two descriptions, the number equations, given by 



on 



Sp 



9/i 



0C 



n 



5n 



3tt 2 



(8) 
(9) 



have to be solved. As an example, the phase diagram is calculated for a fixed density n 
and for a fixed imbalance chemical potential (. To do this, the first number equation (jSJ) 
is solved, given values of (, A and Q, to determine the chemical potential /x. The set of 
values (/!, C, A, Q) is then substituted in the free energy T = Q sp + fin. The minima in 
the free energy landscape determine which state is the ground state of the system, for a 
given imbalance chemical potential ( and a given interaction strength l/kpa s . There are 
three local minima that can be identified in the free energy landscape: the BCS-state (spin- 
balanced superfluid) with A ^ 0, Q = 0, the FFLO-state with A ^ 0, Q ^ and the 
normal state with A = 0. Figure [T] shows that the FFLO-state can indeed be the ground 
state of an imbalanced Fermi gas in three dimensions (at zero temperature). In this figure, 
the free energy of the system is shown, as a function of the bandgap A and the wavevector 
Q, relative to the Fermi energy Ep and the Fermi wavevector kp respectively. 
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FIG. 1: Free energy landscape for an imbalanced Fermi gas in 3D at zero temperature, for different 
values of the imbalance chemical potential: (a) Q/Ep = 0.141 (b) (,/Ep = 0.152 (c) (/Ep = 0.164. 
As the level of imbalance is increased, the competition between the local minima becomes apparent. 
Under the right circumstances, the FFLO-state is the ground state of the system (b). The value 
of the interaction strength used here is tr^ — = — 1. 

For relatively small values of the imbalance chemical potential (, the system is in the BCS 
ground state (a). When ( increases, the system undergoes a first order transition into the 
FFLO-state (b). When the imbalance chemical potential increases further, the system con- 
tinuously goes over into the normal state (c). Figure [T]only shows these phase-transitions for 
one specific value of the interaction strength, near the BCS-limit (1/kpas = —1). However, 
since the first number equation is used to calculate the value of the chemical potential /i, the 
description is also valid for the complete BCS-BEC crossover regime. It must be noted that 
the mean field approximation breaks down in the unitarity limit. However, the FFLO-state 
is expected to form only in the BCS region of the BCS-BEC crossover, where we expect 
the mean field description to give qualitatively correct results. The phase diagram of the 
system is shown in figure [2j This diagram shows that, theoretically, the FFLO-state can be 
formed in an imbalanced Fermi gas in 3D, but since it only occurs on a tiny section of the 
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phase diagram, it may be hard to observe this state experimentally. Our results coincide 



with recent theoretical results 
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FIG. 2: Phase diagram of an imbalanced Fermi gas in 3D at zero temperature, for fixed density 
n. As the imbalance chemical potential £ increases, the system undergoes a first order transition 
from a spin-balanced superfluid (BCS) to the FFLO-state. Above a critical imbalance (which is 
dependent on the interaction strength), the FFLO-state continuously goes over into the normal 
state. 



III. MODELING A ONE-DIMENSIONAL OPTICAL POTENTIAL 

As shown in section [Til the problem with detecting the FFLO-state in an imbalanced 
Fermi gas in 3D is that it exists only in a relatively small section of the BCS-BEC phase 
diagram. In this section, we describe the 3D imbalanced Fermi gas in a ID optical potential. 
There are two main reasons why such a potential can stabilize the FFLO-state. The first 
reason is that in an imbalanced Fermi gas in 3D, the FFLO-state can have a wavevector Q 
in an arbitrary direction. This freedom of choice leads to low-energy bosonic excitations, or 
Goldstone modes, which render the FFLO-state unstable. In the presence of a ID optical 
potential however, it will be energetically favorable for the FFLO-state to form in the direc- 
tion of the optical potential. This will limit the choice for the wavevector Q to just one value, 
thus suppressing the Goldstone modes, which is expected to stabilize the FFLO-state. The 
second reason is that the optical potential will enhance the ID modulation of the FFLO- 
order parameter. We therefore expect the enhancement to be largest when the wavevector of 
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the FFLO-state is equal to the wavevector of the optical potential. The present mean-field 
treatment does not include the effect of excitations such as the Goldstone modes, but it does 
include the energy lowering of the modulated order parameter due to the optical potential. 

In this section we propose two approaches of modeling a 3D imbalanced Fermi gas in 
a ID optical potential. In both approaches, the optical potential is described by using a 
modified dispersion relation. In section [III Al we model the opticaLpotential by introducing 
an anisotropic effective mass in the direction of the potential [34] (from here on this is 
supposed to be the ^-direction). This approximation is valid when the Fermi energy of the 
system lies near the bottom of the lowest Bloch band, i.e. in the case of low density or 



a short-wavelength optical potential. In section IIIIB} we model the optica. 



potential by 



treating the full lowest Bloch band in the tight binding approximation |35l-l38l|. This model 
is valid when the Fermi energy lies in the lowest Bloch band (otherwise more bands have to 
be considered), but contrary to the first case, it does not need to lie at the bottom of the 
band. In both section IHI Al and section [Til Bl the optical potential is supposed not to forbid 
tunneling, as this would inhibit the formation of the FFLO-state. This implies that we will 
treat the imbalanced Fermi gas in 3D in a ID optical potential as a three dimensional system 
with a one-dimensional periodic modulation. 



A. Anisotropic effective mass 

In this section, the ID optical potential is modeled through the use of a modified effective 
mass of the fermionic particles in the direction of the optical potential 

e(k,k z ) = k 2 + ^-, (10) 

where m z is the effective mass of the particles in the ^-direction. Here and for the remainder 
of the paper, k 2 = k 2 + k 2 denotes the magnitude of the in-plane wavevector. This is the 
wavevector which lies perpendicular to the laserbeam which creates the ID optical potential. 
The derivation of the thermodynamic potential for this case is analogous to the derivation 
in section HIl In the present derivation however, it is assumed that the FFLO-state will form 
in the ^-direction Q = (0,0, Q), because this is energetically favorable to other directions, 
due to the anisotropy introduced through (flUj) . The resulting thermodynamic potential, in 
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the limit for temperature going to zero, is given by 
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with the modified notations 



£k = 

CQ,k 



k 2 + 



2m, 



/i 



Q,k 
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(12) 
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■k z Q-C 



The number equations are still given by (jHJ) and (JSJ), but the density n has changed to 



n 



\J2m z 
3tt 2 



(13) 



because of the modified dispersion relation (flOj) . In the limit m z -4-1/2 the thermodynamic 
potential (ITTj) and the density (TT3]) converge to the corresponding expressions in the case 
of an imbalanced Fermi gas in 3D, described in section [Til Expression ([13]) implies that 
when the effective mass m z changes, the density n changes with it. It would be interesting 
however, to compare the phase diagrams for Fermi gasses with different effective masses at 
equal density This can in fact be achieved, because we have found that the thermodynamic 
potential of the system with an effective mass m z can be rescaled to the thermodynamic 
potential of the system with effective mass m z = 1/2 fl6]), using the following scaling relation 

1 \ r-. / .giii 



n 



^u, C, A, Q, m z , ^-j = s/2m z ~Q sp \^t, (, A, 



(14) 



a s J ' X \ /< 2m z 2 a s y/2m z 

From a theoretical point of view, this rescaling property is time-saving for calculations and 
gives a deeper insight into the role of the effective mass m z . The main advantage is however, 
that all physical properties can be studied at the same density. This property relates to 
experiment, because when an external potential is turned on, the effective mass is altered, 
but the average density will remain the same. The effect of changing m z on the BCS-BEC 
crossover phase diagram is shown in figure [31 This figure shows the FFLO phase boundaries 
of the imbalanced Fermi gas in 3D for different values of the effective mass m z , before 
rescaling according to (fl4j) (and hence at different densities). There is a tilting of the FFLO- 
region about a fixed point at unitarity. Furthermore, there is an increase in the width of 
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FIG. 3: Phase diagrams of 3D imbalanced Fermi gasses in a ID optical potential (at temperature 
zero), where the optical potential was modeled by altering the effective mass m z (in the z-direction) 
of the fermionic particles. Here, the phase diagrams are shown for different effective masses m z . 
The shaded regions indicate the FFLO phase. These phase diagrams are shown before rescaling 
according to (|14|) and hence the densities differ for the various systems with different effective 
masses. After rescaling, all FFLO-regions map onto the case of m z = 1/2 (indicated in light gray). 

the FFLO-region (relative to the abscissa) as the effective mass increases. By using the 
scaling- relation (TT4I) . the phase diagrams in figure [3] can be rescaled to equal density. After 
rescaling, the phase diagrams for the different effective masses maps onto the phase diagram 
of the imbalanced Fermi gas with isotropic effective mass [m z = 1/2), described in section 
HT1 From this we conclude that the FFLO-state is not fundamentally influenced by an optical 
potential in which the Fermi energy lies near the bottom of the first Bloch-band. This can 
be explained by the fact that no fundamental anisotropy is introduced into the system by 
altering the effective mass, because independent of the effective mass, the system can be 
scaled back to the case of the imbalanced Fermi gas where the effective mass equals 1/2. 

B. Bloch-dispersion 

In section IIII Al it was shown that a more fundamental anisotropy is needed, in order 
for the optical potential to have an effect on the FFLO-state. In this section, we model 
a ID optical potential in the tight-binding approximation, using the first Bloch-band. For 
this purpose, the quadratic dispersion in the z-direction is replaced by a periodic dispersion 
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k 2 + 5 



1 — cos 



Ql 



(15) 



Here Ql is the wavevector of the optical potential and 5 is a prefactor with units of energy, 
given by [381 ] 



5 



7T 



cxp 



(16) 



with Vo the depth of the potential and Er the recoil energy given by E R = 2ir 2 h 2 /m\ 2 , with 
A the wavelength of the optical potential and m the mass of the fermionic particles. In the 
limit for small k z expression (115D simplifies to ( fTUl) . with 



Ql 

<5tt 2 ' 



(17) 



Given the new dispersion f fl5|) . the thermodynamic potential for this system can be calcu- 
lated. The result is 
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with the following notations 

^Q,k = k 2 + 5 



E* 



1 — COS \ %-q£ \ ( ' ()S 
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I -, ()8 (|iL)co8(^)]-„} + A 2 



(19) 



k 2 + 6 

k C Qlk = C-5sin(|^)sin(^ 

It can easily be shown that expression (fl8l) is equal to the corresponding thermodynamic 
potential ( II ip of the anisotropic effective mass case, in the limit Ql — > oo, <5 — > oo with m z 
held constant, according to (1T7|) . The two number equations again are given by ([HD and ([2]) 
and the density n can be calculated using the general expression 
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dk 



e i 
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(20) 



which yields 
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(21) 
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Here it must be noted that our derivation is only exact if Ep < 25, because otherwise more 
than one Bloch band has to be considered. As in the previous sections, the phase diagram 
for this system can be constructed by studying the local minima of the free energy. Figure H] 
shows a comparison between the phase diagram of an imbalanced Fermi gas in 3D and the 
phase diagram of an imbalanced Fermi gas in 3D subject to a ID optical potential, modeled 
in the tight-binding approximation using the first Bloch band. 
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FIG. 4: Comparison between the phase diagram of an imbalanced Fermi gas in 3D and an imbal- 
anced Fermi gas in 3D subjected to a ID optical potential (with 5 = 0.5 and Ql = 1.2), modeled 
in the tight-binding approximation, using the first Bloch band. Both phase diagrams are at tem- 
perature zero. The stabilizing effect of the optical potential enlarges the FFLO region by a factor 
3 to 6 when compared to the case of the imbalanced Fermi gas without optical potential. When 
C/Ep w 0.48, the width of the FFLO-region is maximal (for this specific choice of 8 and Ql)- This 
resonant enhancement of the FFLO-region occurs when the wavevector of the FFLO-pairs is equal 
to the wavevector of the optical potential. 

This figure shows that the FFLO-region is enlarged by a factor 3 to 6, due to the stabilizing 
effect of the ID optical potential. Figure H] further shows that at C,/Ep ~ 0.48 a transition 
point occurs, where the FFLO-region reaches a maximum width (relative to the abscissa), 
and narrows quickly for larger values of £. This effect finds its origin in the magnitude of the 
wavevector of the FFLO-pairs Qfflo- When the imbalance chemical potential ( increases, 
Qfflo increases likewise to accommodate for the widening gap between the Fermi surfaces of 
the two spin-species. At a certain level of imbalance (in the case of figure |4] at Q/ 'Ep ~ 0.48) 
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Qfflo equals the wavevector of the optical potential Ql- At this point, the FFLO-state 
is optimally enhanced, because the spatial modulation of the FFLO-state is equal to the 
spatial modulation of the optical potential. This results in a maximal width of the FFLO- 
region. When ( increases further, Qfflo retains the constant value Ql, and is not able to 
grow any further. This effect is shown in figure Hence, we can conclude that, although 
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FIG. 5: The wavevector of the FFLO-state Qfflo, as a function of the imbalance chemical 
potential £ (relative to the Fermi energy Ep). The value of Qfflo increases with increasing 
imbalance, until it reaches the value of the wavevector of the optical potential Ql (in this case 
Ql = 1.2), which is the saturation point. The single overshoot point for 5 = 0.6 is probably a 
numerical inaccuracy. 



the imbalance has increased further, the optical potential forces the FFLO-state into a state 
where the form of the FFLO-order parameter matches the form of the optical potential. 
This results in a narrowing of the FFLO-region, because the wavevector of the FFLO-state 
is not sufficiently large anymore to bridge the gap between the Fermi-surfaces of the spin-up 
and spin-down particles. The value of the transition point, where Qfflo becomes equal to 
Ql, rougly increases linearly with the value of 5, as shown in figure [61 Qualitatively, this 
is because the rate of change of the FFLO- wavevector Qfflo with increasing imbalance 
chemical potential (, decreases when 5 becomes larger. This means that ( has to be larger 
(compared to the case of lower 5) for Qfflo to reach the limiting wavevector of the optical 
potential Ql- The advantage of this resonant enhancement of the FFLO-state is that, for 
a given level of imbalance, an optimal stability region for the FFLO-state can be created, 
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simply by tuning the wavelength of the ID optical potential. It should be noted that, when 
considering an imbalance chemical potential ( smaller than 0.48, 5 has to become smaller 
than 0.5 and more bands have to be taken into account in order for our description to be 
exact. We do not treat this case in the present paper. 

To obtain a more direct link with experimental parameter values, we convert our units 
back to SI units for the three situations depicted in figure [6j A possible choice of atoms 
which we considered is m K atoms in a one-dimensional harmonic trap, with a density of 
10 13 cm -3 . For instance, the case with 5 = 0.6 and Ql = 1.2 corresponds to the case of an 
optical potential with wavelength equal to 724 nm, a recoil energy of 309 nK and an optical 
potential depth of Vq/Er ~ 2.22. The numerical values for these experimental parameters 
in the case of 5 = 0.5 and of 5 = 0.7 are depicted in figure M panels (a) and (c) respectively. 
For the illustrative cases of figure [6] we use Ql = 1.2, but theoretically, any choice of Ql was 
possible because we found that the wavevector of the optical potential Ql acts as a scaling 
parameter, according to the following scaling relation: 

Q sp ( /i, C, A, Q, 5, aQ L , j^—) = a Q sp I (, A, -Q, 5, Q L , -j^— J . (22) 

In principle this means that we can vary Ql from zero to infinity. However, there exist some 
limitations on this parameter. First there is a lower limit for Ql because below a certain 
value of Ql no value of the optical potential depth Vo can satisfy equation ( TlBl) . given 
values for 5 and for the recoil energy Er. Second, when the depth of the optical potential 
becomes too large, particles will be confined in the direction of the optical potential, thus 
inhibiting the formation of FFLO-states. This sets an upper limit for the ratio of Vq/Er 
and subsequently for the value of Ql- 

During the course of our work, Loh and Trivedi published their results on the LO-state 
in a 3D cubic lattice 29J]. They found that in a 3D cubic lattice, the LO-state was more 
stable than the FF-state. Since we already find a substantial increase in the FF-state using 
a ID optical potential, we expect that the effect on the LO-state will be similar or larger. 
It would be interesting to apply our ID-potential scheme also to the LO-case. 
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IV. CONCLUSIONS 



We have described the FFLO-state in an imbalanced Fermi gas in 3D within the path- 
integral framework, by choosing a suitable saddle-point at which the atomic pairs have a 
finite centre-of-mass momentum. As a platform to address the case of a 3D imbalanced Fermi 
gas in a ID optical potential and to validate our path-integral description, we rederived the 
zero-temperature phase diagram for an imbalanced Fermi gas in 3D. For this case, our results 
coincide with recent theoretical results. As a proposal to stabilize the FFLO-state we have 
studied an imbalanced 3D Fermi gas in a ID optical potential. This potential was modeled 
in two different ways. For the first model, where we considered anisotropic effective masses, 
we have found that this model is a rescaling of the case of a 3D imbalanced Fermi gas with 
isotropic effective mass. In the second model, we described the effect of the ID optical 
potential using the first Bloch band in the tight-binding approximation. In this case we 
have found a substantial increase in the stability region of the FFLO-state, as compared to 



the case of the 3D Fermi gas without the ID optical 



29 



potential. Related results were recently 
The advantage of our ID optical 



found in the case of a 3D cubic optical lattice 
potential scheme, compared to a 3D cubic optical lattice, is that it allows to find an optimal 
stability configuration for the FFLO-state with a given level of imbalance, by tuning the 
wavelength of the optical potential. This resonant enhancement of the FFLO-region occurs 
when the wavevector of the FFLO-pairs is equal to the wavevector of the optical potential. 
This tunability makes a ID optical potential a suitable experimental configuration for the 
stabilization of the FFLO-state. We therefore propose that this concept can facilitate the 
experimental observation of the FFLO-state in an imbalanced Fermi gas in 3D. 
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l/k F a s 

FIG. 6: Several phase diagrams of an imbalanced superfluid Fermi gas in 3D subjected to a ID 
optical potential, for increasing values of 5. The value of the imbalance chemical potential £ at 
which the maximal FFLO-region occurs, roughly scales linearly with the value of 5. For each 
level of imbalance, an optimal FFLO-region can be found, by tuning the wavelength of the optical 
potential. 
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